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We investigate the photon statistics of light transmitted from a driven optical cavity containing 
one or two atoms interacting with a single mode of the cavity field. We treat arbitrary driving fields 
with emphasis on departure from previous weak field results. In addition effects of dephasing due to 
atomic transit through the cavity mode are included using two different models. We find that both 
models show the nonclassical correlations are quite sensitive to dephasing. The effect of multiple 
atoms on the system dynamics is investigated by placing two atoms in the cavity mode at different 
positions, therefore having different coupling strengths. 
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I. INTRODUCTION 



In this paper we report on extensions to previous work on dynamical cavity QED effects in the photon statistics 
of transmitted light from a driven optical cavity coupled to an ensemble of two level atoms. Much work has been 
done on structural cavity QED effects such as energy level shifts and the modification of spontaneous emission rates. 
These structural effects can be seen to arise from semiclassical models. In addition work has been done on dynamical 
effects where the coupling between the cavity field and atoms has a significant effect on the evolution of the system, 
in particular in the strong coupling regime where a single quantum of energy and hence single quantum fluctuations 
give rise to nontrivial dynamics. In this regime the field can neither be viewed as mildly perturbed by the atoms 
(good cavity limit), nor are the atoms mildly perturbed by the field (bad cavity limit). For a review of the work on 
structural and dynamical effects in cavity QED, see Ref. 

The problem of a single two level atom coupled to a single mode field was originally studied by Jaynes and Cummings 
and extended to many atoms by Tavis and Cummings . These models have been extended in recent theoretical 
work to include spontaneous emission and cavity field decay , and atomic transit time broadening and detunings 
[0. Nonclassical correlations in photon statistics which violate a Schwarz inequality have been predicted for this 
system, including photon antibunching [defined here as g^^^(0)+ > g^^^O)], and sub-Poissonian statistics [g^^^O) < !]■ 
Also other effects have been predicted which we refer to as overshoots and undershoots [|g*-^''(T) — 1| > \g''^\0) — 1| 
where 5'^-'(t) is the normalized second order correlation function]. Examples of these nonclassical correlations from 
previous weak field results are shown in Fig. |l| Figure ^a) shows photon antibunching and sub-Poissonian statistics, 
(b) shows an overshoot violation, and (c) shows an undershoot violation. 

Photon antibunching has been seen experimentally in this system by Rempe et. al. Overshoot violations 
have recently been seen by Mielke, Foster, and Orozco |9j. In general the theory matches the experiments in terms 
of qualitative behavior while the quantitative size of the nonclassical effect does not. This has led us to consider 
complications in the experiments which may be responsible for the discrepancy including deviations from the weak 
field limit, and dephasing due to the atoms entering and leaving the cavity. The experiments use an atomic beam 
to introduce atoms into the cavity so that the time of flight across the mode is on the order of ten spontaneous 
emission lifetimes We expect that dephasing due to atomic traversal of the cavity will have a detrimental effect 
on nonclassical correlations. In addition deviations from the weak driving fleld limit and interactions with 'spectator' 
atoms far from the mode waist may be important. These effects are investigated in this paper by numerically solving 
the master equation for the system and by quantum trajectory simulations. Rather than investigating all possible 
effects at once, we isolate them and try to understand what is most critical. 

The general outline of this paper is as follows. In Section pi w e present the physical model of the system under 
investigation and describe the methods of solution. In Section pi| we discuss the photon statistics of the transmitted 
light outside the weak field limit. Section IV presents two models of atomic transit dephasing and the resulting photon 
statistics. In Section ^ we include effects of a spectator atom with a coupling which is a fraction of the maximum 
coupling, and finally we conclude in Section VI. 
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II. PHYSICAL MODEL 



The system under investigation is an extension of the Jaynes-Cummings Hamihonian which includes effects of 
atomic and cavity field decay as well as a coherent driving field. A schematic diagram of the system is shown in Fig. 
Ij. The field and atomic Hamiltonians are given by 



Hp = huJcCi a (1) 



3 

and the atom-field interaction in the rotating wave approximation is given by 

Haf = ^ i%j f aVi - aa^^j . (3) 



The cavity field creation and annihilation operators are and a respectively and CTj_ and al are Pauli operators for 
the jth two-level atom. The atom-field coupling strength is determined by 

5j = ^ ofc \r si'^ (4) 



\2heoV 

where fj, is the dipole transition matrix element between the two atomic states, V is the cavity mode volume, and 
sin kzj takes into account the position of the atom in the mode. In previous work it was assumed that the atoms sit 
at antinodcs of the field where the coupling is a maximum. Here we allow the atoms to be placed anywhere in the 
mode so that a range of couplings are allowed for different atoms. The cavity field is driven by a classical laser field 
with the driving field-cavity field coupling described by the Hamiltonian 

Hl = iTiE (ate"*""* - ae'""*) (5) 

where E is the classical laser intensity scaled such that E/k is the photon flux injected into the cavity. Throughout 
we assume the atom, cavity field, and driving field are on resonance {loq = uJc = t^a = ^d)- 

Dissipation in the system gives rise to nontrivial irreversible dynamics. Cavity field damping and atomic population 
and polarization decay are described by superoperators acting on the density matrix of the system which are derived 



using standard methods |£0 11 1 . Cavity field damping is described by 

CfP = n (2apa) — a) ap — pa) (6) 



where k is the rate of cavity field damping. Atomic population and polarization decay are described by 

ly 
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where 7 is the spontaneous emission rate of an atom. The full master equation in the Born-Markov approximation is 
then 

p = - J [Ha + Hf + Haf + Hl,p\ + CfR + Cap ^ Lp. (8) 

A numerical solution of the master equation is carried out in the Fock state basis, and also a quantum trajectory 
simulation is developed from the master equation. 

A. Numerical solution of the master equation 

The master equation in the Fock state basis is 

Pn, + -in,+ = -gVn + lpn+l--m,+ - gVm + lpn, + :m+l,- + Ey/npn-l, + -m.+ (9a) 
-^E^/mpn, + -m-l,+ - Ey/jl + lpn+l, + -m,+ - Ey/m + lpn, + -m+l,+ 
+ 2Ky^ {n + l)(m + l)p„+i_ + ;,„+i_+ - [K[n + m) + 7] Pn, + ;7n,+ 
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-'n. — :m,— 



Ti — 1 , + ; m , ■ 



-1,+ + E^/npn-l., 



'Ey/mpn--ni-l- - E^n+ lpn+l-:m,- - E^/rnTlpn~■,m+l. 



(9b) 



Pn,-\-,m,— 



-g\/n+Tp„+i,+;m,- + + -E\/np„-l, + ;m,- 

+-EVmpn, + ;m-l - - -B\/nTTp„+l. + ;,„,- - £;\/m + lpn, + ;m+l, 

+2Ky^ {n + l)(m + l)p„+i,+;m+i,- - [K{n + m) - 7/2] 



3c) 



(9d) 



where Pn,±;m,± = (^1 ilpl'Ti, ±) and + and — denote upper and lower atomic states respectively. 

We have numerically solved the master equation for the steady state for arbitrary driving field by truncating the 
Fock basis at a point where the population of \njnax,i) is less than 10~^. The second order correlation function 



,(2)(^)^ (at(0)at(r)a(r)a(0)) 



(ata)2 



(10) 



is calculated from steady state matrix elements using the quantum regression theorem due to Lax ||l^ . 



B. Quantum trajectory simulation 



We have developed a quantum trajectory simulation of this system from the master equation following the formalism 
of Carmichael |13|. We unravel the master equation into both a piece describing continuous evolution and a set of 
collapse operators in a way which is based on a simulated photon counting experiment. 



Cp= iC-S)p + Sp 



(11) 



where {C — S) p is identified as the terms which can be written as commutators or anticommutators and Sp is identified 
as all terms which can be written as O'^ pO. This particular unraveling is well suited for studies of photon statistics as 
the O's represent quantum jumps due to emission of a photon. The continuous evolution of the system is described by 
[C ~ S) p while Sp describes collapse events which punctuate the evolution. We define a closed system Hamiltonian 
and a dissipative Hamiltonian from the unraveled master equation as 



{C-S)p = --[Hs,p] + [Hd,pI 



(12) 



[Ha + Hf + Haf + Hl,p] 



where O, p denotes the anticommutator of O and p. A non-Hcrmitian Hamiltonian which reproduces the continuous 
L J + 

evolution of the density matrix is defined as 



H = Hs + iTiHo- 



(13) 



The rest of the master equation enters as collapse operators which are applied at random times when i?(0, 1) < Pc 
where R{0, 1) is a random mmiber between zero and one and 



Pc = {^j\6^d\^)dt. 



(14) 



The time step size is {20r)~^ where r is the fastest rate in the problem. In the event that we get two collapse processes 
in a single time step, we use a random number to choose one of the collapses. The time step is chosen so as to minimize 
such occurrences. For this system, the collapse operators are a and ai_ corresponding to emission of a photon from 
the cavity field and the atom respectively. In Section [V, we describe dephasing due to an atom leaving the cavity 
using another collapse operator. 

Because this unraveling of the master equation is based on photon counting experiments, the calculation of the 
second order correlation function is carried out quite naturally. The collapse operator a corresponds to emission and 
detection of a photon from the cavity field mode. We calculate g^^^r) by building up a histogram of delay times 
between photon detections averaged over a long evolution time in a way analogous to experimental measurement. 
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III. NON-WEAK DRIVING FIELD 



The photon statistics of the transmitted Hght have aheady been calculated in the weak field limit using a truncated 
five state basis where the system has up to two quanta of energy in it 1^. The three types of nonclassical behavior 
previously discussed have been seen in subsequent experiments, however these experiments are not strictly in the 
weak field limit. It is of interest then to calculate the photon statistics for arbitrary driving field and to see to what 
extent the nonclassical effects persist. It is expected that for a strong enough driving field, the atoms saturate and the 
nonclassical photon correlations will be washed out because the cavity will basically contain a coherent state which is 
only mildly perturbed by the presence of the atom. 

Looking at the photon correlations in the weak field limit from the point of view of quantum trajectories, we can 
interpret the nonclassical effects as resulting from the collapse of the wavcfunction. The detection of the first photon 
emitted from the steady state collapses the wavefunction of the system (IV'ss) aj-i/'ss)) and the subsequent time 
evolution as the system returns to the steady state determines the photon correlations. The second order correlation 
function is given by the probability of detecting a second photon normalized to the probability of detecting a photon 
in the steady state. Here we start the system in the steady state, collapse the wavefunction, and let it evolve to get 

i 

= (15) 

This assumes that there is usually only one photon emitted as the system returns to steady state which is a good 
approximation in the weak field limit. An example of this is shown in Fig. |^ for the case of the overshoot violation. 
Outside of the weak field limit the photon correlation s are altered for two reasons. Most simply, the time evolution 
following a collapse from steady state will be altered by the stronger driving field. Another effect, however, is the 
presence of multiple collapses before the system returns to the steady state. Consider a multiple collapse process. 
The first photon comes from the steady state and collapses the wavefunction of the system {\ipcoiiapsei) = a|'0ss))- 
Now the time evolution occurs as before. However, the second photon collapses the system to a new state which 
depends on the delay time since the emission of the first photon [\1pcoiiapse2) = o,\ipcoiiapsei{T))]. If a third photon is 
emitted before the system returns to steady state then its delay time will depend on the details of the evolution from 
\ipcoiiapse2)- When averaged over many instances, this process will wash out the nonclassical effects because of the 
different evolution following different possible states, \'ipcoUapse2) ■ An example of this process is shown in Fig. ^ where 
the conditioned photon number undergoes two collapse events. This figure shows three delay times of I/7, 2/7, and 
3/7 with different evolutions of the conditioned photon number resulting in each case. 

Figure H shows the time evolution of {a^a)c following a photon emission from the steady state for a variety of system 
parameters. The overshoot persists in the evolution of the field following emission of a photon from the cavity for 
a driving field as large as E/Esat = 0.41. The undershoot and sub-Poissonian statistics survive for driving fields as 
large as E/Egat — 0.8 and E/Egat — 0.37 respectively. The saturation field strength Egat is the driving field for which 

72 

{n) = Tlsat = (16) 

The photon statistics of the transmitted field are shown in Fig. ^ for the three types of nonclassical effects seen in this 
system at a variety of driving field intensities. Figure ^(a) shows ^'^•'(t) for system parameters (5/7 = 1, K/7 = 0.77) 
which produce an overshoot violation of the Schwarz inequality [g^^Hr) > ^^^•'(O)] in the weak field limit. At a 
driving field of E/Egat = 0.17 the overshoot violation is gone, thus this nonclassical effect is quite dependent on the 
weak driving field. Figure |^(b) shows photon statistics for system parameters {g = 2/7,^/7 = 5) which produce an 
undershoot violation of the Schwarz inequality [1— (7(2) > g(2)(^o) — 1] in the weak field limit. Here the nonclassical 

effect disappears at a driving field of E/Egat = 0.28 showing that this is a more robust effect. Figure ^c) shows 
photon statistics for system parameters (5/7 — 1, k/7 — 1.6) which produce photon antibunching [g^^\0)+ > ^'^•'(O)] 
and sub-Poissonian statistics [g^'^\0) < 1] in the weak field limit. In this case the nonclassical effect persists until 
E/Egat = 0.16 where the system shows slight bunching and super-Poissonian statistics. (Notice that the nonclassical 
effects are not as robust as the time evolution of the cavity field would indicate. Therefore the destruction of 
nonclassical effects are in part a result of multiple photon processes.) For all system parameters the transmitted light 
becomes super-Poissonian as the driving field is increased. 



4 



IV. ATOMIC TRANSIT DEPHASING 



We now turn our attention to the effects of atomic traversal of the cavity on the photon statistics. In previous work 
it has been assumed that the atoms are all fixed at antinodes of the cavity field and so had the maximum coupling 

1/2 

5o = M (^c/27ieoy) . Experiments on this system have used atomic beams to send atoms through a cavity. This 
atomic traversal of the cavity will introduce two new effects. First, the atom/cavity field coupling will depend on the 
position of the atom in the cavity. One might think that this would destroy the nonclassical correlations. However, 
the atoms with the largest coupling interact most strongly with the field and are most likely to contribute to the 
correlations. So the atoms near an antinode will have the largest contribution and other atoms may have little effect 
on the correlations. This issue will be further addressed in Sec. The second effect of atomic traversal is dephasing 
which occurs when an atom enters or leaves the cavity. It is this effect which we consider in this section. 

We have used two approaches to model the dephasing due to atomic traversal. The first is to add a term to the 
master equation which describes nonradiative decay of atomic polarization. 

P = Cp + 7p/i {(TzPCTz - p) (17) 

This term in the master equation has its origins in coUisional processes and so may or may not accurately describe 
the dephasing which occurs when an atom leaves the cavity. 

The second approach uses a quantum trajectory simulation of the system to model the dephasing. In this approach 
we assume that there is always exactly one atom in the cavity. An atom leaves the cavity and another atom enters the 
cavity in the ground state at a rate ^ph- We can assume the atom enters the cavity in the ground state, but it is not 
immediately clear how to deal with the state of the exiting atom. This atom is in some superposition of excited and 
ground state and these states are entangled with the cavity field state. One approach would be to leave the photon 
number distribution of the cavity field unchanged using a collapse operator which has the following action on the 
state of the system: 

IV-) =5](ce,n|e,7i)+cg,„|g,n))^ ItA,) =^(c2„+4„)|5,n). (18) 

n n 

However, this is not a consistent application of the quantum trajectories. Consider the evolution of the atom after it 
leaves the cavity. The atom at some later time may emit a photon into the vacuum meaning it was in the excited state 
when it left the cavity; or it will never emit a photon, meaning it was in the ground state when it left the cavity. In 
general, the atom and environment and by entanglement the atom/cavity system will then be described by a density 
operator. However, we wish to use a pure state to describe the atom/cavity system conditioned on the detection of 
transmitted photons. To be consistent we must use a pure state to describe the atom after it has left the cavity. We 
use a collapse operator which picks either the excited state field distribution or the ground state field distribution of 
the system and then places the new atom in the ground state. This operator has the following action: 

IV') = X! ('^£'"1'^' + Cg,„|5, n)) I^Z-c) = ^ Ce,„|5, n) with probability c^ .^ (19a) 

n n 

= ^ (ce,„|e,7i) +Cg.n\g,n)) IV'c) = '^Cg,n\g,n) with probability Cg„. (19b) 

n n 

This collapse operator is then applied at a Gaussian distributed series of times with average \/^ph and a full width 
of as this approximates the traversal times of atoms with a Maxwell-Boltzmann velocity distribution. This 

model of dephasing differs from coUisional dephasing in two important ways. First, it does not enter the deterministic 
Hamiltonian evolution between collapses at all, whereas the coUisional dephasing in the trajectory picture would have 
both a collapse component as well as the decay of coherence in the continuous evolution. Second, this dephasing 
always places the atom into the ground state whereas the coUisional dephasing places the atom in the ground or 
excited state with probabilities determined by the relative populations. 

Now we turn to our results and a comparison of the two types of dephasing. Because the traversal dephasing does 
not affect the deterministic evolution of the system we expect that for a given dephasing rate it will be less destructive 
of the nonclassical photon statistics than coUisional dephasing. In Fig. |7| we show the second order correlation 
function with coUisional dephasing. All three types of nonclassical effects are quite sensitive to this dephasing with 
7p/i = 0.05 destroying the overshoot [Fig. 0(a)] and undershoot [Fig. |^(b)] while the sub-Poissonian statistics survive 
until ^ph = 0.2 in Fig. ^(c). Our results for the transit dephasing are shown in Fig. ^. The overshoot in Fig. ||(a) is 
again extremely sensitive to dephasing with classical statistics for ^ph/l — 0-05 while the sub-Poissonian statistics in 
Fig. |8|(b) are more robust, surviving up to 'fpuH — 
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V. TWO ATOM EFFECTS 



We now consider the effect of placing two atoms inside the cavity either both at antinodes of the field, or allowing 
one of the atoms to be arbitrarily placed so that its coupling is in the range to (70 . Dephasing is not considered in 
this Section. In the experiments conducted on this system it is likely that there is some effect from 'spectator atoms' 
which are located away from an antinode of the field and so do not contribute to the nonclassical photon statistics. If 
there are enough of these atoms they may simply absorb light and emit it out of the cavity thus effectively decreasing 
the quality of the cavity. As a first step towards understanding the effect of spectator atoms, we place one extra atom 
in the cavity with a coupling which is some fraction of the original atom's coupling. We use a quantum trajectory 
simulation to calculate photon statistics. 

In Fig. H we plot ^'■^''(t) while allowing the coupling of the spectator atom (172) to vary from 50/IO to go- We see 
that for all three sets of parameters the photon statistics vary continuously from the single atom result to the two 
atom result with no qualitative deviation in the photon statistics. This suggests that the spectator atoms can have an 
observable effect on the statistics, but they do not destroy the nonclassical correlations at least when close to the ideal 
condition of having a single atom in the cavity at a time. A group of many spectator atoms may be more detrimental 
to nonclassical correlations. 

We now consider the case in which both atoms are placed at antinodes of the field. The photon statistics for this 
system have been solved for an arbitrary number of maximally coupled atoms using a set of symmetrized Dicke states 
to describe the atomic excitation |7|. This would correspond to an experimental setup where it is not possible to 
tell which atom spontaneously emitted a given photon. However, there would be experimental situations where the 
symmetrized states are not valid states with respect to spontaneous emission. Here, using quantum trajectories, we 
consider both a symmetrized collapse and an unsymmetrized collapse for spontaneous emission. The operator for the 
Dicke collapse is 



so that when a spontaneous emission event occurs, both atoms are collapsed symmetrically. For the non-Dicke collapse 
the atomic operators are used separately so that one atom or the other collapses. 

For weak driving field we expect that there will be no difference between these types of collapse because in this limit 
we only detect photons emitted from steady state. However, for stronger driving fields we begin to detect photons 
emitted from the collapsed state and the two collapses give a different collapsed state {a^\ipss) or a-tlipss) versus 
Cpickel'fpss)), and therefore, different photon statistics. 

In Fig. |l^ g^'^\T) is plotted for the two types of collapse for a driving field E — 0.5. Figure ^(a) shows a 
significant difference as the nonclassical statistics are completely gone for the non-symmetrized collapse. Figure p^(b) 
shows no dependence on the type of collapse. Figure ^(c) shows a mild dependence on the type of collapse with a 
slightly larger value of g^^^ (0) for the non-symmetrized collapse. 



We have investigated extensions to previous theoretical work on a driven atoms/cavity system with dissipation. 
We have calculated the normalized second order correlation function for the transmitted light including effects of 
arbitrary driving fields, non-radiative dephasing, and arbitrary coupling strength for multiple atoms. We have found 
that nonclassical field states are easily destroyed by deviations from the weak field limit and by non-radiative dephasing 
modeled as both coUisional dephasing and as atomic transit dephasing. We have also found that allowing two atoms 
in the cavity with different atom/field coupling strengths does not have a detrimental effect on the nonclassical field. 
The experiments which have been done on this system have not really in the weak field limit. As i? — *■ the number 
of counts also goes to zero so it is difficult to carry out experiments in this regime, but this work suggests that it is 
important. 



The authors would like to thank Howard J. Carmichael, Luis Orozco, and Greg Foster for many useful conversations. 
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FIG. 1. Examples of nonclassical photon statistics in the weaJi field limit, (a) shows sub-Poissonian statistics and photon 
antibunching. (b) shows an overshoot violation of the Schwarz inequality, (c) shows an undershoot violation of the Schwarz 
inequality. 

FIG. 2. A diagram of the system. The cavity field decays at a rate n, the atom spontaneously emits at a rate 7, and there 
is an electric dipole coupling between the atom and cavity field with a strength g. 

FIG. 3. The conditioned photon number normalized by the steady state photon number for the cavity field. This is identical 
to 5^^^(r) in the weak field limit. Parameters are g/7 = 1, /t/7 = 0.77, £/7 = 0.01. 

FIG. 4. The conditioned photon number normalized by the steady state photon number for the cavity field. In these plots, 
two collapses occur giving rise to different evolutions depending on the delay between collapses. The delays between collapses 
are (a) I/7, (b) 2/7, and (c) 3/7. On the average, processes such as these will wipe out the nonclassical correlations. 

FIG. 5. Time evolution of the conditioned cavity photon number following emission of a photon from the cavity and collapse 
of the wavefunction from steady state. The plots are for (a) gj-y = 1, K/7 = 0.77, -B/7 = 0.1 (solid line), 0.2 (dashed line), 0.3 
(dotted lino), (b) g/j = 2, k/j = 5, E/-y = 0.1 (solid line), 0.5 (dashed line), 1 (dotted line), (c) 5/7 = 1, /t/7 = 1.6, E/'y = 0.1 
(solid line), 0.5 (dashed line), 1 (dotted line). 

FIG. 6. Photon statistics of the transmitted field for varying driving field strength. The plots are for (a) 
5/7 = 1,k/7 = 0.77, E/j = 0.025 (solid line), 0.125 (dashed line), 0.2 (dotted line), 0.35 (dash-dotted line). (b) 
g/7 = 2,fv/7 = 5, -E/7 = 0.025 (solid line), 0.25 (dashed line), 0.35 (dotted line), 0.5 (dash-dotted line), 1 (small-dashed 
line), (c) 5/7 = 1,k/7 = 1.6, S/7 = 0.025 (solid line), 0.25 (dashed line), 0.425 (dotted line), 0.6 (dash-dotted line). 

FIG. 7. Photon statistics of the transmitted field with collisional dephasing. The plots are for (a) 

g/j = 1, k/7 = 0.77,5/7 = 0.1,7ph/7 = (solid line), 0.05 (dashed line), 0.2 (dotted line). (b) 

5/7 = 2,k;/7 = 5,E/y = 0.1,7^^/7 = (solid line), 0.05 (dashed line), 0.2 (dotted line). (c) 
5/7=1, /t/7 = 1.6, £;/7 = 0.1, jph/j = (solid line), 0.1 (dashed line), 0.2 (dotted line). 

FIG. 8. Photon statistics of the transmitted field with atomic transit dephasing. The plots are for (a) 
5/7 = l,/t/7 = 0.77, S/7 = 0.1,7ph/7 = 0.05 (solid line), 0.1 (dashed line), 0.5 (dotted line). (b) 
g/j = l,/t/7 = 1.6,£;/7 = 0.1,7p/i/7 = 0.05 (solid line), 0.1 (dashed line), 0.5 (dotted line). 
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FIG. 9. Photon statistics of the transmitted field with two atoms in the cavity at arbitrary coupling strength. 
The plots are for (a) c?i/7 = 1,k/7 = 0.77, S/7 = 0.1,32/7 = 0.1 (solid line), 0.5 (dashed line), 1 (dotted 
line). (b) pi/7 = 2,^/7 = 0.77, £;/7 = 0.1,52/7 = 0.2 (solid line), 1 (dashed line), 2 (dotted line). (c) 
51/7 = 1,k/7 = 1.6, E/ J = 0.1,32/7 = 0.1 (solid line), 0.5 (dashed line), 1 (dotted line). 



FIG. 10. Photon statistics for the transmitted field with two atoms comparing the symmetrized Dicke spontaneous emis- 
sion with non-symmetrized spontaneous emission. The plots are for (a) 171/7 = 52/7 ~ l,i^/"t = 0.77, £'/7 = 0.1. (b) 
51/7 = 52/7 = 2,k/7 = 5,E/y = 0.1. (c) 51/7 = 32/7 = 1,k/7 = 1.6, E/^ = 0.1. The solid line is the non-symmetrized 
collapse and the dashed line is the symmetrized collapse. 
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